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ABSTRACT 

We present a detailed analysis of the spatially and spectrally resolved 12 C0 J=2— 1 
and J=3— 2 emission lines from the TW Hya circumstellar disk, based on science verifi- 
cation data from the Atacama Large Millimeter/Submillimeter Array (ALMA). These 
lines exhibit substantial emission in their high-velocity wings (with projected velocities 
out to 2.1km s _1 , corresponding to intrinsic orbital velocities >20km s _1 ) that trace 
molecular gas as close as 2AU from the central star. However, we are not able to 
reproduce the intensity of these wings and the general spatio-kinematic pattern of the 
lines with simple models for the disk structure and kinematics. Using three-dimensional 
non-local thermodynamic equilibrium molecular excitation and radiative transfer calcu- 
lations, we construct some alternative models that successfully account for these features 
by modifying either (1) the temperature structure of the inner disk (inside the dust- 
depleted disk cavity; r < 4AU); (2) the intrinsic (Keplerian) disk velocity field; or (3) 
the distribution of disk inclination angles (a warp). The latter approach is particularly 
compelling because a representative warped disk model qualitatively reproduces the ob- 
served azimuthal modulation of optical light scattered off the disk surface. In any model 
scenario, the ALMA data clearly require a substantial molecular gas reservoir located 
inside the region where dust optical depths are known to be substantially diminished in 
the TW Hya disk, in agreement with previous studies based on infrared spectroscopy. 
The results from these updated model prescriptions are discussed in terms of their po- 
tential physical origins, which might include dynamical perturbations from a low-mass 
companion with an orbital separation of a few AU. 
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Introduction 



In the standard model for star formation, the collapse of a slowly rotating molecular cloud 
core at least partially conserves its ang ular momentum by formin g ; a compact, flattene d disk that 
channels mass onto a central protostar (ICassen fc Moosmanlll98ll : iTerebey et al.lll984l ). Once the 
remnant core material is accreted or dispersed, the gravitational potential of the star determines the 
kinematic properties of its disk. At that time, the intrinsic disk velocity field should be described 
well by a Keplerian pattern of differential rotation in circular orbits, v^(r) = Vk = yjGM/r, where 
(f) is the azimuthal direction in a polar coordinate system, G is the gravitational constant, r is the 
distance from the star, and M is the mass enclosed within r. In most cases, the disk-to-star mass 
ratio is low enough that M « M* is a reasonable approximation. Imaging spectroscopy data taken 
with millimeter interferometers first verified these signatures of disk rotation using the optically 



thick emission lines of the carbon monoxide (CO) molec ule (IKoerner et al.1 Il993l : iDutrey et al 



1994 : IKoerner fc Sargentl Il995l : Mannings fc Sargentl 119971 ). Subsequent work demonstrated that 
the rotation patterns were inde ed Keplerian (v^ oc r~°- 5 ) and could be used to make dynamica l 
estimates of central star masses (jGuilloteau fc Dutreylll998l : lDutrey et al.lll998l : ISimon et al.ll2000l ). 



In practice, spectral imaging observations measure the intrinsic velocity field projected onto 
the sky, v \^ s = i^sini, where i is the angle between the disk rotation axis and the observed line 
of sight (e.g., i = 0° corresponds to face-on). A Keplerian velocity field can be exploited to probe 
gas a t small radii in the disk, even if the angular resolution of the observations is relatively limited 
(e.g., llDutrey et al.ll2008l ). This kind of physically motivated super-resolution can be achieved with 
a focus on the intensities (and ideally morphologies) in the high-velocity wings of emission lines. 
For optically thick lines like those produced by the rotational transiti ons of CO, the line strengt h 
depends on the product of the gas temperature and emitting area (e.g.. iBeckwith fc Sargentl 1 19931 ). 
Since the projected emitting areas are typically small in the inner disk, emission in the line wings is 
generally weak; often well below the sensitivity thresholds for most millimeter- wave interferometers. 

To further complicate issues, there are theoretical mechanisms and disk properties that could 
modify this simple prescription for interpreting the projected disk velocity field. For example, in a 
very massive disk, gravitational instabilities could excite spiral waves that drive significant (radial) 
streaming motions. This was suggest ed as a po tential explanation for the complex kinematic 
environment of the AB Aur disk flLin et all 12009 ). alt hough others have associated its (po ssibly) 
non-Keplerian motions with envelope contamination (IPietu et al J 12005 ; ICorder et al.ll2005l ). Gas 
pressure gradients also lead to deviations from Keplerian orbits: the s ub-Keplerian shift exp ected 



from thermal (hydrostatic) pressure is too small to observe directly (jWeidenschillind Il977), but 



magn etic pressure might have a more substantial influence over gas motions (e.g. Jshu et al 



20071, 



I2OO8I ) . Alternatively, a disk might appear to have non-Keplerian motions (even if = Vk exactly) 
if its line-of-sight orientation (i) varies with radius - that is, if the disk structure is warped. 

Given that emission line wings are typically expected to be weak, disentangling these subtle 
(hypothetical) modifications to a simple projected velocity field from the intrinsic properties of the 
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inner disk makes for a difficult task in data analysis. Nevertheless, the potential for this unique 
access to the spatio-kinematic properties and physical conditions of the gas in the innermost regions 
of protoplanetary disks is extraordinarily compelling. Perhaps the best opp ortunity to explore 



these features is with the massive, gas-rich disk around the nearest (d « 54 pc; Ivan Leeuwenl 120071 ) 
classical T Tauri star, TW Hya. Aside from the enhanced sensitivity and spatial resolution afforded 
by its proximity, th e TW Hya disk co i ncident ally ha s the added benefits of a nearly face-on viewing 



geometry (i w 6-7°; lKrist et alJl200(j ; iQi et all 12004 : [Hughes et all 1201 This means that the line 



wing emission generated by gas in the inner regions of the TW Hya disk has relatively small observed 
Doppler shifts from the systemic velocity. Because the low disk inclination angle minimizes radial 
projection effects, subtle departures from a simple (unwarped, Keplerian) model for the projected 
velocity field should be more easily recognizable observationally. 

In this article, we take advantage of the dramatically improved sensitivity now available in the 
Science Verification (SV) data products from the Atacama Large Millimeter/Submillimeter Array 
(ALMA) to help characterize the spatio-kinematic morphologies of the spatially and spectrally 
resolved CO J— 2— 1 and J=3— 2 emission lines from the TW Hya circumstellar disk. A brief 
overview of the data and its calibration are presented in §2. The data is investigated in detail in 
§3, using emission line radiative transfer calculations and considering various toy models for the 
disk properties. The results are synthesized in the larger context of the structure and kinematics 
in the inner regions of the TW Hya disk in §4, and summarized in §5. 



2. Data and Calibration 

The TW Hya disk was observed with nine 12-m elements of ALMA as part of the commissioning 
and SV effort during construction. These data were obtained in two (contiguous) ^4. 5-hour obser- 
vations, using receiver bands 6 and 7 on 2011 April 20 and 22, respectively, with antenna spacings 
providing baseline lengths of ^20-100 m (one antenna was not available for band 7 observations at 
this time). The ALMA correlator was configured to simultaneously observe four 0.5 GHz- wide spec- 
tral windows (two in each sideband) with a 122 kHz channel spacing in each window. One of those 
windows was centered near the 12 CO J— 2— 1 line (230.538 GHz) in band 6 and the 12 CO J=3— 2 
line (345.796 GHz) in band 7, providing a velocity spacing of 0.16 and 0.11km s -1 , respectively. 
The ALMA correlator Hanning smoothes the data, making the effective spectral resolution coarser 
(^0.3 and 0.2km s _1 in bands 6 and 7). Two other spectral windows were designed to provide 
^1 GHz of continuum bandwidth in each receiver band; the remaining windows contained weaker 
spectral lines that are not of interest here. Observations cycled between TW Hya and the nearby 
complex gain calibrator J1037-295 every ^10 minutes. The bright quasar 3C 279 was observed as 
a bandpass calibrator, and brief observations were made of Titan for absolute flux calibration. 

These data were calibrated in the CASA software package (v3.3), carefully following the detailed 
processing scripts kindly provided by the ALMA science verification team. Since those instructions, 
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along with fully calibrated measurement sets, are publicly available onlinq^l, we do not repeat the 
details here. The continuum-subtracted, calibrated spectral visibilities for each CO emission line 
were Fourier inverted (assuming natural weighting) with modest spectral averaging, deconvolved 
with the CLEAN algorithm using a simple polygon mask applied to all channels, and restored with 
a synthesized beam of dimensions 2"8 x 2"4 (at a P.A. of 44°) for J=2-l and 1"7 x 1"5 (P.A. 
= 22°) for J=3— 2. We utilized two different channel spacings: (1) "native", meaning 0.20 and 
0.12 km s _1 for J— 2— 1 and J=3— 2, respectively, to highlight the key kinematic features (note 
that this over-samples the true spectral resolution; this is the averaging advocated in the ALMA 
calibration tutorials), and (2) "binned", meaning 0.35km s _1 for both lines, to simplify comparisons 
and sample models closer to the true Hanning-smoothed velocity resolution. The resulting spectral 
image cubes are indistinguishable from the reference datasets provided by the science verification 
team, with one small distinction: we elected to force the reduced spectral datasets to have a channel 
precisely centered on the known systemic LSR velocity of TW Hya, +2.87 km s _1 . 

The CO J— 2— 1 and J=3— 2 line emission from the TW Hya disk are shown as channel maps 
in Figures [1] and El respectively, with a logarithmic color-scale stretch to emphasize emission over 
a wide dynamic range (a scale bar is marked on the top of each figure). Measurements of the RMS 
noise levels, peak and integrated line intensities for each transition are compiled in Table [TJ The 
channel maps clearly show a resolved pattern of rotation, from the blueshifted northwest to the 
redshifted southeast (with a m ajor axis position an gle of 151°), around the systemic velocity and 



centered at the stel lar position (IQi et al.ll2004l. 120061). W hile these are not the highest spectral (see 



Hughes et al.ll201ll ) or spatial (see lAndrews et al.ll2012l ) resolution views of CO emission from the 



TW Hya disk, they are by far the most sensitive (^10 x higher S/N). This enhanced sensitivity 
reveals a surprising spectral feature: significant (>5 a) emission in the wings of both CO transitions 
out to projected velocities at least ±2.1 km s _1 from the line centers. 

We will demonstrate in §3 that these line wings are primarily generated well inside a radius 
of ^5AU in the TW Hya disk and, coupled with the emission morphology near the systemic 
velocity, they are difficult to reproduce for a disk in Keplerian rotation with a nearly face-on 
viewing geometry and smoothly-varying temperature profile. Because of the value of this relatively 
faint emission at high velocities in what follows, we conducted a variety of tests to ensure that these 
line wings were not artificially generated in the calibration process. First, we verified that the disk 
emission is not spectrally coincident with CO features in the Earth's atmosphere. Weak and narrow 
CO absorption (^lMHz wide; ^10-20% depths) was noted and nagged in the observations of the 
bandpass calibrator 3C 279: however, they are located >20MHz from the TW Hya line center, and 
are therefore well outside the spectral range of interest here. A second test focused on the remote 
possibility that the application of the frequency-dependent system temperatures (T sys ) could have 
artificially broadened the lines. In these data, a coarse measure of the spectral T sys behavior was 
determined directly from observations of TW Hya. The concern is that the bright CO emission 



2 https : // almas cience .nrao . edu/alma-data/science-verif ication . 
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Fig. 1. — Channel maps of the CO J— 2— 1 line emission from the TW Hya disk, with a "native" 
velocity-spacing of 0.20 km s _1 . A logarithmic color-scale bar (top) indicates the line intensities, 
with a vertical line representing the 5 a level (see TableHJ. Each panel is 12" on a side, corresponding 
to a projected physical scale of ^650 AU. Cross-hatches mark the stellar position and orientation 
of the disk major axis (151° east of north). The synthesized beam dimensions are drawn in the 
bottom left panel; LSR velocities (in km s _1 ) are marked in the bottom right of each panel. 

from the disk might produce a weak, but broad feature in the T sys profile that would then be 
propagated into the data as faint emission wings. To test this, we re-calibrated the data instead 
with a spectrally flat, band-averaged T sys correction. This had no discernable effect on the results; 
the observed line wings persist at the same intensity levels. Finally, we investigated whether the 
continuum subtraction process could be responsible. A wide range of techniques were attempted, 
varying the polynomial order and channel ranges employed to determine the underlying continuum 
level on either side of the emission line. In all cases, we found no significant difference in the 
resulting line wings (indeed, they are even present without any continuum subtraction applied). To 
summarize, we find no obvious artificial mechanism in the data acquisition or calibration procedures 
that can produce these substantial emission line wings with the same velocity- widths in two CO 
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Fig. 2. — Same as Figure [H but for the CO J=3— 2 line. 

transitions observed at different times and with different receivers and spectral setups. 



3. Analysis 

Having presented the ALMA data (see Figures [1] and [2]) , we now turn to its interpretation 
with a special emphasis on the resolved spatio-kinematic patterns of the CO emission lines. First, 
we briefly describe a generic framework for producing models of the CO lines in the TW Hya disk 
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(§3.1). Next, we demonstrate that standard models are unable to reproduce the ALMA data with 
the typical assumptions about disk temperature profiles and the spatial variation of the disk velocity 
field projected on the sky (§3.2). Then, we present three options for simple modifications to those 
assumptions that are able to reconcile the data and models (§3.3): a large temperature increase 
inside the TW Hya disk cavity (r < 4 AU; §3.3.1); a deviation from a Keplerian velocity field that 
produces relatively enhanced rotation rates at small disk radii (§3.3.2); or a warp that produces a 
relatively higher line-of-sight inclination angle in the inner disk (§3.3.3). We focus on a qualitative 
illustration of how the data reveal novel spatio-kinematic information on gas at small radii in the 
TW Hya disk, rather than employing a quantitative model-optimization (detailed fitting) approach: 
the latter will be more useful in the near future, when ALMA provides TW Hya data with further 
improvements in dynamic range and angular resolution (as well as sensitivity). 



3.1. Modeling Overview 



We defined a simple, azimuthal- and mirror-symmetric parameterized model for the disk gas 
densities on a cylindrical grid with coordinates (r, 0, z), such that 



p(r, z) 



2ttH 7 



exp 



(i) 



where E and H p are the radial surface density and pressure scale height profil es, respectively. For the 
form e r, we adopted the simil arity solution for a simple viscous accretion disk (jLynden-Bell fc Pringle 



1974 ; lHartmann et al.lll998l . see also Andrews et al. 2009, 2010), 

2-7' 



E(r) = E c - 
r c 



exp 



(2) 



where 7 is the gradient parameter that describes the power-law radial variation of the disk viscosity, 
r c is a characteristic scaling radius, and S c = e • E(r c ). Since we considered only emission from the 
CO molecule, the models were parametrically defined in terms of the CO column density, 



iVco(r) 



jiniH 



E(r), 



(3) 



where X co is the CO abundance fraction in the gas, (i = 2.37 is the mean molecular weight of the 
gas, and trh is the mass of a hydrogen atom. The scale height profile was set by the balance of 
gravity and thermal pressure in the disk, such that 



k b r 3 T(r) 



1/2 



(4) 



where c s is the local sound speed, Q is the Keplerian angular velocity, and is the Boltzmann 
constant. In this formulation, we implicitly assumed that the disk is vertically isothermal (at each 



-8- 



radius), with a radial temperature profile T(r) = T 10 (r /10 AU)~ g , where T 10 = T(10 AU). The disk 
kinematics are determined by a velocity field in the model plane, v = (vr^v^^Vz). In the models 
considered here, v r — v z — 0, and the azimuthal component of v is described by Keplerian rotation, 



v (t>( r ) = v k = y—^- 3 ( 5 ) 

assuming implicitly that the disk mass is negligible compared to the stellar mass. 

The structure model described in Equations (l)-(5) is fully characterized by a set of 6 param- 
eters, {7V C , r c , 7, Tio, q, An additional 3 parameters are required to convert this structure 
to a synthetic CO emission line model projected into the sky plane: a turbulent velocity dispersion 
(£), and a disk viewing geometry specified by the inclination (i) and major axis position angle (PA). 
For a given set of 9 parameters, we used the three-dimensional non-local t hermodynamic equilib- 



rium (non-LTE) molecular excitation and radiativ e transfer code LIME (vl.02: lBrinch fc Hogerheijde 



20 id ) and the data from the LAMBDA database (jSchoier et al.ll2005l ) to calculate level populations 



and generate high-resolution synthetic spectral cubes for the CO J— 2— 1 and J=3— 2 emission lines. 
Those cubes were Fourier transformed and re-sampled onto the spectral and spatial frequency grids 
employed in the ALMA observations. The resulting model spectral visibilities were then passed 
through the same CASA imaging processes that were performed on the data (see §2). 



3.2. The Failure of Standard Models 



We first constructed a "fiducial" model using the standard parameterization describe d in §3.1, 



motivated by the results of our previous analyses of CO em ission from the TW Hya disk flQi et al 



2004 120061 : iHughes et al.ll2008L l201ll : I Andrews et al.ll2012l ). Because of the time-consuming nature 
of the LIME calculations for any given model, we adopted a simplistic exploration of parameter-space 
that utilized the modified Levenberg-Marquardt minimization algorithm MPFIT (jMarkwardtl 120091 ) 
to help find parameters that produced a low x 2 in the spectral visibilities for both CO transitions. 
In that process, we forced the basic parameters that describe the disk densities {7V C , 7, r c } and 
turbulent velocity width {£} to be the same for both CO transitions, but allowed different values 
for the temperature profile parameters (T in, q} for the two emission lines. Th e fiducial stellar 



mass was fixed to M* = 0.8 M (e.g., IWichmann et al.lll998l ; IHughes et all 1201 ll ). and the major 



axis position angle was set to 151° (east of north); the disk inclination {i} was considered a free 
parameter (albeit the same for both CO transitions). Guided by the minimization results, the 
model parameters were slightly adjusted manually to best reproduce the observed channel maps 
(at native resolution). The fiducial model parameters we derived ar e provided in Table p (not 
surprisingly, they are similar to those for the Model sA constructed by lAndrews et al.ll2012l ). The 
model channel maps are compared directly with the data in Figure [3l imaged on the binned velocity 
grid with the same channel-spacing of 0.35 km s _1 for both CO transitions, for clarity. 



2 N C is defined as the CO column density at the characteristic radius, r c ; see Equation (3). 
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Fig. 3. — Comparisons of the observed CO J— 2— 1 and J=3— 2 channel maps (top rows) with the 
fiducial model and the three other standard models described in §3.2 (see Table [2] for parameter 
values). The observed and model channel maps were synthesized on the same binned 0.35km s _1 
velocity gridding for ease of comparison. The map sizes, intensity scale bars, LSR velocity labels, 
and synthesized beam dimensions are the same as in Figures [1] and El 



The fiducial model is a good match to the ALMA data near the systemic velocity, but notably 
fails to reproduce the emission line wings that are observed for both CO transitions. Those wings 
are observed at velocities at 2.1km s _1 (and beyond, albeit at lower statistical significance) from 
the line core, corresponding to intrinsic disk velocities in excess of 20 km s _1 considering the inferred 
disk inclination of ^6°. For this stellar mass and assuming Keplerian orbits, those velocities are only 
reached at disk radii <2AU. Since these CO lines have such high optical depths, their intensities 
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scale with the product of the gas temperature and the projected emitting area (IBeckwith fc Sargent 



19931 : iDutrey et al.l 119941 ). With such a small area available to produce emission at these line wing 
velocities, this discrepancy between the fiducial model and the data is perhaps not so surprising. 
Broadly speaking, there are two approaches that could help reproduce the observed CO emission 
line wings: the first deals with adjusting the physical conditions of the model at small radii, and 
the second relies on modifying the model parameters that control the projected velocity field. 

For the first (structure) approach, a representative model needs to incorporate higher gas 
temperatures at small disk radii (because the optical depth is so high, adjustments to the CO 
column densities have little effect on the line wing emission strengths). Sticking with the standard 
prescription for modeling this kind of data (described in §3.1), that requirement effectively amounts 
to adopting a steeper temperature profile. We constructed such a "high-g" model with a larger 
radial temperature gradient (q = 0.65); parameter values are listed in Table[2j and synthetic channel 
maps are compared directly with the data and fiducial model in Figure [31 By increasing the gas 
temperatures inside ^5AU by roughly a factor of 3, this kind of model is able to reproduce the 
observed intensities in the emission line wings. However, the required tempertaure profile is too 
steep to properly account for the emission morphology generated at relatively large disk radii, and 
therefore ends up being a very poor match for the CO emission near the systemic velocity. 

The alternative (kinematic) approach relies on increasing the projected area of the high- velocity 
disk region responsible for producing the CO line wings. Put simply, the observations demand an 
increase in the projected velocity field (^ bs — v^sini). This can be achieved either by scaling 
up the intrinsic disk velocities (v^ = Vk oc y/M*) and/or the projected line of sight (sini). To 
demonstrate that, we constructed two additional models using the standard prescription outlined 
in §3.1. First, a "high-M*" model employed a larger stellar mass (M* = 1.5 M ) to speed up the 
Keplerian orbital velocities (ity), thereby increasing the extent of the high- velocity emitting area. 
And second, a similar effect was achieved for a "high-z" model that adopts a more inclined viewing 
geometry (i = 8°). Both of these models increase the projected line of sight velocities by a factor 
of ^1.4, which corresponds to a factor of ^3.5-4 increase in the emitting area that produces the 
high- velocity line wings. The parameters for these models are listed in Table El Their corresponding 
synthetic channel maps are compared with the data and fiducial model in Figure El As expected, 
the high-M* and high-i models successfully reproduce the observed CO line wings: however, they 
consequently do a very poor job accounting for the morphology near the systemic velocity. 

In practice, the high-g, high-M*, and high-i models that follow the standard formalism outlined 
in §3.1 all fail to reproduce the observed CO emission line morphologies in the same way^ No 
reasonable combinations of g, M*, and i can simultaneously reproduce both the CO line wings and 
the emission morphology near the systemic velocity. Uniformly scaling up the temperature profile 



3 From a statistical perspective, these kinds of model are formally much poorer matches to the data than the 
fiducial model: their discrepancies with the data are primarily in the spectral channels that contain most of the 
emission signal, and therefore most of the weight (contribution to % 2 ) m any quantitative assessment of fit quality. 



gradient or projected velocity field at all disk radii cannot explain the patterns of the ALMA data. 



3.3. Alternative Models 

The standard modeling formalism explored in §3.2 is unable to explain the observed spatio- 
kinematic emission morphologies of the CO emission lines generated by the TW Hya disk. However, 
that modeling effort demonstrated clearly the "sense" of this failure. To produce the emission line 
wings, we need a model with either much higher temperatures or increased projected velocities 
in the inner disk. However, the observed emission morphology near the systemic velocity was 
explained best with a fiducial temperature structure and projected velocity field. This implies that 
we require a deviation from the standard structure and kinematics model that varies with radius 
in the disk. As before, these modifications to T(r) or ^ bs( r ) — v^smi can be manifested in a 
deviation to the form of the temperature profile, the intrinsic disk velocity field, ity(r), and/or the 
disk viewing angle, z(r), so long as those functional forms decrease with r. Three such alternatives 
were explored here, using simple parametric adjustments to the standard model prescription. 



3.3.1. A Hot Inner Disk 



The first approach we adopted for these alternative models was to implement a structural 
modification to the radial temperature profile. As was shown in §3.2, models with a hot inner disk 
(steeper gradient q) are able to reproduce the observed CO line wing intensities, but the outer disk 
needs to have a more shallow temperature profile to account for the emission morphology near the 
line core. In principle, these features could be reconciled by permitting the gradient parameter q 
to vary with radius. However, such a modification is not so intuitive; instead, we elected to simply 
permit a multiplicative scaling of the temperatures inside a given radius, such that 

T(r) = (ST) ■ T w (r/10AUH (6) 

where ST is a free (constant) parameter in the inner disk, but ST = 1 when r > r^ ot . In effect, this 
model adjustment preserves the fiducial outer disk temperature profile, but artificially scales up the 
temperatures at small disk radii to produce the CO line wings. Based on this modified framework, 
we developed a representative model for the CO emission following the procedure described in §3.2. 
The resulting parameters for this "hot" model are listed in Table El Synthetic observations of the 
model are compared with the data as channel maps in Figure [H as well as in position- velocity 
diagrams (constructed from a slice along the disk major axis) in Figure El 

We find a good match to the CO data by modifying the fiducial model described in §3.2 with 
a scaling factor ST — 3 inside a radius r^t — 4AU. This choice for r^t was n °t arbitrary; it was 
selected because t he TW Hya disk i s known to have substantially diminis hed dust optical depths 
inside that radius (jCalvet et al.ll2002 ; Hughes et al. 2007 : Arnold et al. 2012 ). Notably, the implicit 
assumption for this scenario is that a substantial gas reservoir is still present inside the dust cavity. 
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Fig. 4. — A channel map comparison of the data, fiducial model, and the alternative models 
described in §3.3 (similar to Figure [3|). Model parameters are listed in Table O 



3.3.2. A Non-Keplerian Velocity Field 

The second alternative we considered was a kinematic modification to the intrinsic velocity 
field of the gas in the inner disk. Again, as demonstrated in §3.2, the outer disk velocities can not 
also be scaled up without producing large deviations in the observed emission morphology near the 
systemic velocity. To accomodate both requirements, we introduced a simple scaling factor for the 
disk velocity field, v^(r) = fv^ (see Equation 5). The scaling factor / describes the model deviation 
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from Keplerian orbits and is permitted to vary with radius like a power-law in the inner disk, 



where {r^, x} are free parameters that describe a break radius and the gradient of the scaling 
function, respectively. In the standard models described in §3.1 and 3.2, / = 1 at all disk radii 
(effectively, x = in those cases). Adding in this parametric non-Keplerian scaling factor as a 
slight modification to the standard modeling framework, we constructed representative models for 
the CO emission lines using the procedure outlined above. The resulting parameter values are listed 
in Table EJ the model is compared with the ALMA data in Figures H] and El 

Adopting r\) « 60 AU and x = 0.15 to describe /, we were able to reproduce well the observed 
CO line wings without sacrificing the fit quality of the emission morphology near the systemic 
velocity. The super-Keplerian rotation in the inner disk means that the CO emitting area at the 
observed velocities of the line wings is substantially enhanced relative to the fiducial model (by a 
factor of ~5), and therefore able to generate high- velocity emission at roughly the observed levels. 
The {r^, x} values adopted here are not unique (nor optimized); a smaller break radius and larger 
gradient (or vice versa) could probably accomodate the data with a similar fit quality. However, 
these values and this kind of model modification are sufficient to demonstrate the key point: a 
non-Keplerian deviation in the disk velocity field which decreases with radius successfully accounts 
for the observed spatio-kinematic CO structure emitted by the TW Hya disk. 



We explored another kinematic alternative model that incorporates a warp in the vertical 
structure of the gas disk. Rather than raise the intrinsic orbital velocity of the gas (ity), this model 
increases the observed velocity by changing the line of sight (sini) as a function of radius. In this 
way, a warp - or a disk inclination that increases toward the central star - raises the projected 
velocities in the inner disk, effectively increasing the emitting area at a given observed velocity. If 
the inclination is a decreasing function of r, this model can be designed to have minimal impact 
on the kinematic pattern for the outer disk, promoting a favorable emission morphology near the 
systemic velocity. To construct a warped disk model, we introduced a power-law profile for the 
observed inclination angles, 



aligned with the disk major axis (assuming a vertically thin disk; see the Appendix), where {io, y} 
are parameters that describe the inclination normalization (at 7™o) and its radial gradient, respec- 
tively. For a given set of structure parameters, the model densities and temperatures are rotated to 
this modified viewing geometry. The velocity field is tilted using a quaternion transformation, as 
described in more detail in the Appendix. Incorporating this parametric formulation for a warped 




(7) 



3.3.3. A Disk Warp 




(8) 
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Fig. 5. — CO position-velocity diagrams constructed from a cut along the disk major axis (P.A. = 
151°) for the data and the fiducial, hot, non-Keplerian, and warped models described in §3.3. For 
reference to gauge the extent of the line wings, the 5 a brightness level is marked with a contour. 

disk viewing geometry, we constructed models for the TW Hya CO emission lines as before. The 
resulting parameter values for a representative warped model are included in Table El synthetic 
observations of the models are compared with the data in Figures [4] and El 

Models that include a modest warp, with iq « 8° at ro = 5AU and a gradient y = 0.15, 
provide an excellent match with the ALMA data for both CO emission lines, reproducing well the 
observed line wings and emission morphology near the systemic velocity. Like the non-Keplerian 
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model described in §3.3.2, this warp formulation successfully produces high-velocity line emission 
by increasing the projected CO emitting area at the observed wing velocities by a factor of ^5 (by 
design). And also like that case, the warp model parameters we determine or even the formulation 
we adopt are not unique, nor optimized. Instead, they are intended to conceptually verify that a 
warped disk geometry can reproduce in detail the fundamental features of the observed CO line 
emission with the addition of only a single model parameter; an inclination gradient, y. 



3.4. Additional Comments 



It is worth pointing out that the alternatives discussed in the previous subsections are in no 
way mutually exclusive: various combinations of these effects could be used to achieve similar 
results. Moreover, there are additional alternatives to a simple disk model that could presumably 
account for these ALMA observations, but have not been explored in detail here. Perhaps most 
compelling among these is a model that explicitly incorporates a spatial variation in the broadening 
term of the line profile. For example, it seems reasonable that a model with a turbulent line width £ 
that decreases with radius should be able to re-distribute some emission into the high-velocity line 
wings obser ved here. The £ values we adopted are commensurate with the turbulence constraints 
imposed by iHughes et al.1 (120111 ) for the outer disk, made using observations with ^4-5 x higher 
spectral resolution than the ALMA data. However, we cannot rule out higher £ values at small 
disk radii. We attempted a simple modification to our fiducial model that permits t he turbulent 
velocity widths to scale with the sound speed of the gas, using the lHughes et al.l (120111 ) upper limit 
to normalize £(r) in the outer disk: but given the low £ at large r and modest radial variation of c s , 
this has little effect on the CO line wings. If a variable turbulent line width parameter is the origin 
of these wings, £ must increase more steeply than c s in the inner disk. Given the limited velocity 
resolution of these ALMA data relative to the subsonic turbulent velocity dispersions in the TW 
Hya disk, we felt it was premature to push the further exploration of such models. Nevertheless, it 
remains an interesting possibility that could be investigated with future observations. 



4. Discussion 



We have conducted a detailed analysis of the 12 CO J— 2— 1 and J=3— 2 line emission from the 
gas-rich disk around the nearby young star TW Hya, making extensive use of the science verification 
data from the Atacama Large Millimeter Array (ALMA) commissioning effort. With the substantial 
improvement in line sensitivity afforded by ALMA (even during its construction), we identified a 
novel feature in these CO spectra: the presence of relatively faint line wing emission extending out 
to high projected velocities of at least ±2.1 km s _1 from the (systemic) line center, roug hly twice 



the maximum velocities previously inferred for the T W Hya disk from less sensitive data (lOi et al 



2004, 


2006; 


Hughes et al. 


2011: 



of the TW Hya disk, gas on Keplerian orbits that generates these line wings must be located within 
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a few AU from the central star. This kinematic super-resolution was employed in an exploration 
of the temperature and velocity structure of the inner disk, utilizing a non-local thermodynamic 
equilibrium molecular excitation and radiative transfer code and a standard modeling prescription 
for circumstellar disks. We found that the typical assumptions of a simple, Keplerian disk with 
a smooth gas temperature profile were unable to simultaneously produce both the CO line wings 
and the observed line emission morphology near the systemic velocity. Instead, we developed three 
alternative model possibilities that are able to qualitatively account for the ALMA data by making 
small (parametric) adjustments to the spatio-kinematic structure of the gas in the inner disk. 

In one such alternative, we considered a parametric sc aling of the gas temperatures inside the 
known dust-depleted cavity at the disk center (r < 4 AU; ICalvet et al.l 12002 : iHughes et al.1 120071 ) 
in an effort to boost the CO line wing intensities without substantially disturbing the emission 
morphology near the line core. With the outer disk structure based on a fiducial model, an increase 
of the gas temperatures in the dust cavity by a factor of ^3 successfully accounts for the ALMA 
CO observations. While such a temperature scaling might seem simplistic and artificial, it is in fact 
consistent with continuum radiative transfer models for the TW Hya disk: the low optical depths 
inside th e disk cavity mean that stellar radiatio n is effectively unattenuated, leading to substantial 
heating 



Andrews et al.ll2012l ). If we assume that gas and dust temperatures are 
correlated (if not identic al), then the factor ST « 3 derived here is quite reasonable. For reference, 



the models developed by I Andrews et al. have dust temperatures that increase by a factor of 

^4 over a narrow (^1 AU-wide) region around the cavity radius (see their Figure 3). It is interesting 
to consider the key implication of this model: the CO emission line wings detected by ALMA may 
offer an indirect (and not very practical) observational signature of a low-density dust cavity, but 
they also require that this cavity still hosts a significant reservoir of molecular gas. 

Of course, the presence of gas inside the TW Hy a disk cavity has been inferred from other di- 
agnostics: accretion indicator s (IMuzerolle et al.ll2000f). as well as emiss ion lines of simple atoms and 
molecules (INajita et al.ll2010l ). including CO (e.g., IRettig et al.ll2004l ). Based on a rotational dia- 
gram constructed f rom t he infrared spectrum of CO fundamental (v = 1—0) rovibrational emission 
lines, ISalyk et al.1 (120071 ) used a simple slab model to estimate T = 650-1500 K and N co = 0.2- 
1.4 x 10 19 cm -2 for an emitting region that corresponds to r « 0.1-1 AU. Taking the mean tem- 
perature from the two CO lines, we find T = 640-1660 K and N co = 14-140 x 10 19 cm -2 for our 
"hot" model in the same region - remarkably similar, keeping in mind that column density es- 
timates from optically thick lines are inherently ambiguous. However, other analyses of the CO 



rovibrational spectrum suggest much low er temperatures (by a factor of ^2; IRettig et al.l 120041 ; 
Pontoppidan et al.ll2008l : ISalyk et al.ll2009l ). Moreover, our "hot" model assumes enhanced temper- 
atures over a continuous emitting area out to 4 AU: if a smaller area were adopted, the temperature 
scaling factor 5T would have to be proportionately increased to compensate for the observed line 
wing intensities (notably to uncomfortable levels that start to impinge on the thermal dissocia- 
tion energy). Given the uncertainties involved in analyzing both kinds of CO data, it is not clear 
whether a "hot" model derived from the rotational lines (as in §3.3.1) is necessarily in conflict 
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(i.e., too hot) with the rovibrational spectrum. In any case, it is interesting to note that these two 
complementary methods could eventually be used to obtain independent constraints on the gas 
properties inside the TW Hya dust disk cavity. In the near future, ALMA can provide sensitive 
high angular resolution images of these same CO transitions that will resolve the spatio-kinematic 
morphology of their line wings. That information will be crucial for locating the CO emitting area, 
enabling more direct comparisons between the infrared and millimeter-wave diagnostics. 

While the first model described above was devoted to modifications of the structure in the 
inner disk, the remaining two alternatives focused on adjustments to the disk kinematics. Rather 
than increasing the intrinsic line intensities from a fixed emitting area, these other models produced 
the observed CO line wings by increasing the effective emitting area itself. In one such scenario, 
we modified the intrinsic velocity field in the inner disk with a simple parametric deviation from 
Keplerian orbits. We found in §3.3.2 that a model which incorporated super-Keplerian rotation 
in the inner disk successfully reproduced the observed CO line emission. Although there is con- 
siderable uncertainty involved, our "non-Keplerian" model rotated ^30% faster than Keplerian 
speeds at r = 10 AU, and ^80% faster at 1 AU (/ « 1.3 and 1.8, respectively). If this deviation 
was restricted to much smaller disk radii only, the corresponding scaling factor / would need to 
be substantially larger to account for the observed CO line wings. Admittedly, of the three al- 
ternative models explored here, this scenario is the least physically-motivated. While in principle 
a (reverse) pressure gradient could accelerate the dis k material, most models that consider suc h 
gradients instead call for sub-Keplerian motions (e.g.. Iweidenschilling 1977 : Shu et al. 2007 . 2008 ). 
However, it is possible that non-azimuthal gas velocities could mimic the parameterization used 
here. For example, a molecular outflow with high mass-loss rates, directed primarily in the vertical 
direction (a v z component) and launched from the inner disk, mig ht play a role. The TW Hya 
disk is thought to drive a substantial photoevaporative wind (e.g., Pascucci et al.l [20T1I ) . but the 
molecular properties of that flow have not yet been explored in sufficient detail to easily compare 
with observations. Another example might be a veloci ty field with radia l streaming motion (a v r 
component), perhaps directed along spiral arms (e.g., buillen et aDbood ). No such structural fea- 
tures have yet been identified for the TW Hya disk; but there are similar disks that do exhibit this 



20031 : iFukagawa et al.1 12004 12006 



kind of asymmetry ([W einbe rger et al.lll999l:lAugereau et al.lll999l : iGrady et al.ll200ll ; IClampin et al 



Muto et al 



201^ K Though these possibilities seem exotic, they 



remain in consideration until realistic, testable models for more complex disk velocity fields are 
developed. 

The final alternative model we developed to account for the ALMA CO data also relied on 
modifying the observed disk velocity field, v \^ s = v^sini. In this scenario, we maintained the 
standard assumption of Keplerian orbits (v^ = v^) and instead permitted the projection factor 
sini to vary with radius. In §3.3.3, we developed a simple parametric description of a disk warp 
that simultaneously accounts for the CO line wings and the nearly face-on emission morphology 
near the systemic velocity. That simple model has an inclination of ^4° in the outer disk and 
increases only slightly to ^10° at r = 1 AU; the mean line-of-sight viewing angle averaged over the 
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disk area is comparable to the ^4 -7° inferred previously fcrist et alJboOO ; Weinberger et al. 2002 : 
Qi et al.ll2004l : lHughes et al.ll201ll ). Warps have be en directly imaged in scattered light for the edge- 
on (e ffectively gas-free) debris d isks around /3 Pic fjKalas fc Jewittlll995l : LHeap et al.ll2000l ) and AU 
Mic (jLiull2004l ; iKrist et al.ll2005l ), and indirectly inferred from quasi-periodi c photometric variations 
in op tical /infrared light curves in othe r edge-on cases, notably KH 15D (Ichiang fc Murray- Clay 
2004 ) and AA Tau (IBouvier et al.l ll999). Warped viewing geometries have also been speculated for 



a few additional cases, in cluding the disks around GM Aur (from CO spe ctral imaging; iDutrey et al. 



1998 : iHughes et all 12003 ) and HD 100546 (from scattered light images: iQuillenl 12009 ). 



In the specific case of the TW Hya disk, iRoberge et al.l (120051 ) have cited potential evidence for 



a warp based on an azimuthal asymmetry in their wideband Hubble Space Telescope (HST)/STIS 
coronagraphic images of optical light reflected off the disk surface. They identify a roughly si- 
nusoidal variation of scattered light with the disk position angle in a region ~ 1.3-1. 6" from the 
central star, with a peak near P.A. = 234° and a trou gh near 54° (m e asured E of N, and notably 



roughly ±90° from the disk major axis orientation). Roberge et al. f 20051 ) speculated that this 
azimuthal modulation in the brightness profile might be the consequence of a varying illumination 
pattern in the outer disk, generated by the shadowing effect of a warped disk geometry at much 
smaller radii. This is effectively the same type of model we have employed to explain the ALMA 
observations of the CO fines, s o nat urally we were motivated to test this possibility against the 



behavior of the IRoberge et al.1 (120051 ) scattered light data. To do so, we implemented a simplified 



dust structure based on the "warp" model described in §3.3.3, assuming that the dust and gas are 



Calvet et al. 


2002; 


Hughes et al. 


2007) 



re lative mass rati o , grai n composition, and size distribution used in the disk atmosphere models 



of I Andrews et al.l (120121 ). The three-dimensional Monte Carlo radiative transfer code RADMC-3! 
(v0.30: C. P. Dullemond) was then used to generate a hig h-resolution synthetic image at a wave- 
length of 0.8 /im, assuming isotropic scattering. Following IRoberge et al.l (120051 ) . we calculated a 
normalized intensity profile in 20° azimuthal bins in an annular ring from r = 70-88 AU. 

The results are shown in Figure [6] and demonstrat e a striking resembla nce to the azimuthal 
modulation pattern inferred from the HS T/STIS data (IRoberge et all 120051 . see their Figure 9). 
Taken at face value, this remarkable agreement might be considered compelling independent ev- 
idence for a very modest warp in the TW Hya disk structure. Again, this scattered light model 
is not a fit to the HS T/STIS observations of scattering off dust grains, but rather the predicted 
behavior based on one model that is able to describe ALMA observations of the CO gas structure. 
The sinusoidal pattern of the azimuthal profile is a result of material in the inner disk - where the 
warp is enhanced - shadowing dust at larger radii. The phase of that pattern is set by the assumed 
warp (rotation) axis, which was arbitrarily aligned with the observed major axis of the disk pro- 
jected on the sky (151° E of N). However, in practice the phase depends only relatively weakly on 



4 http : //www . ita . uni-heidelberg . de/$\sim$dullemond/sof tware/radmc-3d/ 
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Fig. 6. — A comparison of the (normalized) azimuthal profile of optica l light scattered off t he TW 
Hya disk surface from both HST/STIS coronagraph data (red circles; iRoberge et al.ll2005l ) and a 
prediction based on our "warp" model (black curve; note that this is not a fit to the scattered 
light data). Both profiles were derived by averaging the intensities in 20° bins from an annulus 
extending from r = 70-88 AU. The sinusoidal modulation is a natural consequence of shadowing 
by a warped structure in the inner disk (higher inclination angles at small radii), which produces 
a varying illumination pattern of the disk surface at larger radii. 



the warp axis orientation: shifting the latter by ±20° leads to only modest changes to the pattern 
phase (i.e., the data are not quite able to quantify the warp orientation). The pattern amplitude 
depends more sensitively on the adopted model parameters, as it is set primarily by the height of 
the optically thick dust layer in the inner disk. If the inner disk were substantially cooler, and 
therefore had lower scale heights, the level of shadowing would decrease and the amplitude of the 
azimuthal scattering pattern would be diminished. Indeed, if we instead adopt the cooler, midplane 
dust temperatures found by ([Andrews et al.ll2012l ) to set the dust scale heights using Equation 4, 
we find a modest (^20%) decrease in the peak amplitudes that are more consistent with the scat- 
tered light data. Perhaps more important in this case is the truncation of the dust distribution for 
r < 4 AU. If we were to assume no dust depletion in the disk cavity, the shadowing of the outer disk 
would be much more pronounced, with a normalized amplitude of ^1 in the azimuthal variation 
of scattered light. In essence, there is a substantial degeneracy between the spatial distribution of 
dust in the inner disk (i.e., the cavity size and depletion factor, as well as the local scale heights) 
and the parameters that describe the warp. Nevertheless, for a set of simple ass umptions grounded 
in the ALMA CO data (and the dust cavity from the SED and resolved images; ICalvet et al.ll2002l : 
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feature of the iRoberge et al 



Hughes et al.ll2007l:lAndrews et al.l 2012 ). we find a warp model is remarkably consistent with a key 



scattered light data. 



Warps have been p ostulated to ori ginate from a wide range of phenomena, including g ravita- 
tional instabilities (e.g., Sellwoodl 2010). intense radiation field s ( Armitage fc Prinde 1997 ). star- 



disk magnetic interactions (e.g., iTerquem fc Papaloizou 



dynamical perturbations from a stellar flyby However 



2000 



Flahe rty fc 



Muzerolle 



and 



none 



of these seem particularly plausible in this case, given the known properties of TW Hya and its disk. 
In principle, a massive planet embedded in the inner dis k with an orbit inclined out of the disk 



plane could also generate the modest warp inferred here (ILubow fc Ogilvidl200ll : iMarzari fc N elson 
2009 ) . While this sc enario has been extensively explored for the f3 Pic debris disk (jMouillet et^l 



19971 : iDawson et all 1201 ll ). the high gas densities in the TW Hya disk could lead to the substantial 
damping o f a planet inclination - a nd therefore the warp structure - on a relatively short (viscous) 
timescale ([Papaloizou fc Lin 1995 ). However, given the large uncertainties on the disk viscosities 
and densities, as well as the range of potential companion properties (e.g., masses, orbits), there 
is still a reasonable likelihood that a companion could maintain the modest orbital inclination re - 
quired to sustain a small warp over the lifetime of the TW Hya disk (e.g., iBitsch fc Kleyl l201lh . 
Moreover, other effects related to planet-p lanet scattering in the syst em could produce repeated 
dynamical excitations of the gas disk (e.g., iThommes fc Lissauerll2003l ). 



W e would be remiss not to point out that there is an alternative explanation for the lRoberge et al 
(hood ) scattered light asymmetry: namely, a spatial variation in the scattering phase function of 
the grains. This latter possibility might better account for the observed kinks in the radial scat- 
tered light profiles, which incidentally our toy model does not reproduce well (although perhaps 
could, given a considerable modeling effort beyond the scope of this article). Other studies of the 
TW Hya scattered light disk from HST m ake no clear mention of an azimuthal asym metry (e.g.. 



Krist et al.1 l2000l : IWeinberger et al.ll2002l ). The original WFPC2 discovery paper bv iKrist et al 
fcood ) does note that the "right" (we assume west) side of the disk appears brighter (consistent 
with the Roberge et al. asymmetry), but they suggest that the discrepancy co uld just as easily be 



explained by an instrumental or calibration artifact. We should also note that IPontoppidan et al 



(120081 ) used infrared spectro-astrometry measurements of the CO fundamental rovibrational lines 
to determine an inclination of only ^4-5° in the inner ^0.1-0.5 AU of the TW Hya disk. The data 



prese nted here are not directly sensitive to such small radial scales, but the IPontoppidan et al 
(j2008l ) claim of a warped outer disk (where i is an increasing function of r) is clearly not consistent 
with the ALMA data. These inferred variations in inclination all could be compatible if the disk 
tilt is not monotonic with radius, but instead has some maximum inclination at a few AU. In any 
case, a simple warp model is a compelling potential explanation for both the ALMA CO spectral 
images and the HST/STIS scattered light data, but it is not necessarily the correct (or unique) 
solution. 



In principle, all three of the alternative models that successfully explain the ALMA CO data 
could ultimately be produced by dynamical interactions between the disk and a companion located 
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a few AU from the central star. Nevertheless, despite the known dust depletion at small radii or 
potential evidence for a warp in this disk, there is not yet any firm evidence for a close companion 
to TW Hya. The most definitive high-resol ution contrast lim its for this system rule out companion 
masses >O.O14M at 2-4 AU separations fjEvans et al.ll2012l ). suggesting that any perturber would 



need to be substellar or pla netary in nature. A previous claim fo r a "hot Jupiter" (jSetiawan et al 



20081 ) has been refuted (IHuelamo et al.ll2008l : lRucinski et al.ll2008l ). and, although there is a sugges- 



tive asymmetry in new mid-infrared data that is consi stent with a faint po int source near the edge of 
the disk cavity, its interpretation remains uncertain (jArnold et al.ll2012l ). Along with deeper com- 
panion searches in the near future, we expect that some of the techniques presented here could be 
used to help provide new, empirical constraints on how dynamical interactions affect the kinematic 
and structural properties of the TW Hya gas disk and others like it. 



5. Summary 



We have studied the spatio-kinematic properties of the CO gas in the disk around the nearby 
young star TW Hya, making use of the ALMA science verification datasets provided during the 
facility construction. While these spectral images of the CO J— 2— 1 and J=3— 2 emission lines 
from ALMA have only modest angular and spectral resolutions (and dynamic range), they provide 
roughly an order of magnitude improvement in sensitivity wh en compared to previous datasets 
fbi et all 120041 . l2006l : IHughes et all boill : Undrews et~aD I2OI2I ). The key results of our analysis 
include: 



1. We find clear evidence for substantial CO emission in high- velocity line wings, extending out 
to at least ±2.1 km s _1 from the systemic velocity. Assuming Keplerian rotation for a nearly 
face-on disk viewing geometry, the projected velocities of these wings correspond to intrinsic 
disk velocities >20km s _1 , indicating that they trace molecular gas at disk radii < 2-5 AU. 

2. Simple models that follow the standard prescription for analysing optically thick molecular 
lines are unable to simultaneously reproduce both the CO line wings and the spatially resolved 
emission morphology near the line cores. Instead, the ALMA data suggest that the inner 
regions of the TW Hya disk must have either much higher gas temperatures and/or a faster 
projected velocity field compared to what would be expected from emission in the outer disk. 

3. Guided by those findings, we developed three alternative models that reproduce the ALMA 
data well by including small parametric modifications to the standard modeling prescription. 
One makes a substantial adjustment to the disk structure by scaling up the gas temperatures 
inside the known dust-depleted cavity in the TW Hya disk (r < 4AU). Two others include 
alterations to parameters describing the observed disk velocity field: the first by permitting 
the intrinsic velocities in the inner disk to increase to super-Keplerian rates, and the second 
by allowing the line-of-sight viewing angle of the disk to increase close to the central star. 
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The latter model of a disk warp finds some independent support in its ability to account for 
the azimuthal asymmetry noted in an optical scattered light image (IRoberge et al.ll2005l ). 



4. Although no one model is explicitly favored by the current suite of ALMA data, they all 
include one implicit requirement: that there be a significant reservoir of molecular gas inside 
the ^4 AU-radius region in the TW Hya disk known to have diminished dust optical depths. 



Perhaps most important, this study illustrates the impending capability of ALMA for providing 
a substantial advance in our understanding of molecular gas in protoplanetary disks. In principle, 
the powerful sensitivity of ALMA and the kind of kinematic super-resolution highlighted here will 
offer novel access to the characteristics of gaseous material in the innermost regions of these disks. 
In the near future, those capabilities will be supplemented with additional detail, as ALMA pushes 
toward significant improvements in both dynamic range and angular resolution as well. 



We are grateful to Rebekah Dawson, Diego Muhoz, Christian Brinch, and Joanna Brown 
for their helpful suggestions, as well as the ALMA commissioning and science verification team 
for kindly providing detailed instructions that guided the calibration of these data. A.M.H. is 
supported by a fellowship from the Miller Institute for Basic Research in Science. The Atacama 
Large Millimeter/submillimeter Array (ALMA), an international astronomy facility, is a partnership 
of Europe, North America and East Asia in cooperation with the Republic of Chile. This article 
makes use of the following ALMA Science Verification data: ADS/ JAO.ALMA#2011. 0.00001. SV. 



A. Notes on Three-Dimensional Rotations with Quaternions 

In the general case, the rotation of a three-dimensional vector quantity can require some un- 
wieldy algebraic manipulation. A much more elegant and straightforward technique to accomplish 
the same result relies on the use of the quaternion algebra developed by W. R. Hamilton in the 
1840s! Although quaternions have been used in many aspects of physics (notably in the Dirac 
formulation of quantum mechanics) and are commonplace in fields like aerospace engineering and 
computer graphics, we have only rarely come across them in the astrophysics literature. If the warp 
structure postulated in §3.3.3 turns out to be a common disk feature, it might be useful to provide 
a brief tutorial on the use of quaternions in three-dimensional vector rotations. 

Quaternions are a generalization of complex numbers. Any quaternion q is represented by the 
combination of a scalar s and a vector v = (v x , v y ,v z ), such that 

q= [ s ,v\ = [s,v x ,v y ,v z ]. (Al) 



5 The original series of articles - entitled On Quaternions, Or On a New System of Imaginaries in Algebra and 
published in installments from 1844-1850 - is a recommended resource; it has been transcribed and edited by David 
R. Wilkins and graciously posted online: see http://www.emis.ams.org/classics/Hamilton/OnQuat.pdf 
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Quaternion algebra has non-commutative multiplication 

P * Q — [ r 5 u] * [ 5 5 ^] — [ r s — u • v, rv + su + uxv], (A2) 

an inner product, p • q = [r, u] • [5, v] = r s + u-v, and conjugation, q = [s, — v] = [s, — i^, —v y , —v z ]. 
For the specific task of applying a rotation, the definitions of an angle 9 and a rotation axis a can 
be compactly encoded in a quaternion, where 

q= [cos (0/2), a sin (0/2)]. (A3) 

Any vector v can be properly mapped into its rotated counterpart w by generating corresponding 
quaternions with zero scalar components, vo = [0, fT] and wo = [0,1?], and then applying a simple 
quaternion algebra manipulation, 

wo = q * vo * q. ( A4) 



A Worked Example: To demonstrate more concretely how the formalism described above can 
be used in practice, we describe the methodology used in §3.3.3 to apply a warp rotation to the 
model disk described in §3.1. For simplicity, we assume that the rotation axis corresponds with the 
Cartesian x-axis, so that a = (1, 0, 0) by definition. In our models, this axis is by default aligned with 
the observed major axis position angle, oriented 151° E of N projected on the sky. We considered 
a rotation by an angle i(r), which depended on the distance from the model center, r — y / x 2 + y 2 
(see §3.3.3; Equation 8). In that context, consider a parcel of gas located at v = (x p ,y p ,z p ) where 
x p — z p — {) and y p = 1 AU (i.e., 1 AU along the y-axis, perpendicular to the warp rotation axis, 
and in the disk midplane). The quaternion vq = [0,1/] describes the parcel position. Likewise, the 
warp - a radially varying rotation - is also encoded in a quaternion, q = [cos (z/2), sin (i/2), 0, 0]. 

The warp can be applied following Equation A4, and will result in a rotated position for the 
initial gas parcel that corresponds to the vector component of the quaternion 

w = [cos (i/2), sin (i/2), 0, 0] ★ [0, 0, y p , 0] ★ [cos (i/2), - sin (i/2), 0, 0] 
= [0,0,y p (cos 2 (i/2) - sin 2 (i/2)), 2y p cos (i/2) sin (i/2)] 

= [0,0,y p cosi,y p sini], (A5) 

where the final step employed elementary trigonometric identities to recover the simple result 
expected for a cell originally located orthogonal to the rotation axis. Using the parameters in 
Table [2] to reconstruct the angle i, the updated location of the gas parcel after rotation is described 
by wq = [0,0,0.98 AU, 0.18 AU]; that is, slightly closer to the rotation axis and now lying at some 
modest height (z = 0.18 AU) above the midplane. Repeating this approach for all model locations 
initially in the midplane (z p = 0) produces the warped structure sketched in Figure [3 

Of course, the algebra involved in reproducing Equation A5 is substantially more complex for a 
more general location if, validating the simplicity in the quaternion approach encapsulated in Equa- 
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Fig. 7. — Schematic of a flat model grid (left) in the disk plane (z — 0) and its warped counterpart 
after a quaternion rotation (right). Given the modest rotation postulated in §3.3.3, the grids are 
shown with an extreme aspect ratio of 1:1:20 to highlight the subtle warp about the x-axis. 

tion A4. For interested users, there are established and well-documented quaternion manipulation 
packages available online for both IDlJ§ and pythoE@ (among other common coding packages). 



e http: //cow. physics .wise . edu/$\sim$craigm/idl/math.html 



'http : //exc . harvard . edu/mta/ASPECT/tool$_-$doc/pydocs/Quaternion . html 
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Table 1. Synthesized Channel Map Properties 



line 




noise (1 a) 


peak l v 






[km s -1 ] 


[Jy beam -1 ] 


[Jy beam -1 ] 


[Jy km s -1 ] 


12 C0 J =2-1 


0.20 


0.009 


7.8 ± 0.8 


17.5 ±1.8 




0.35 


0.008 


7.3 ± 0.7 


17.5 ±1.8 


12 C0 J=3-2 


0.12 


0.020 


9.7 ±1.0 


37.2 ±3.7 




0.35 


0.015 


10.6 ± 1.1 


36.0 ±3.6 



Note. — Properties of the synthesized channel maps at native and 
binned spectral resolutions. The noise is defined as the RMS variation 
per channel in line-free regions. The peak I v refers to the peak intensity 
in a channel of fixed width Av. The integrated flux densities {F v — 
J Ijydv) correspond to the sum of the emission at all velocities in a 
5" x 5" square box centered on the stellar position. The quoted errors 
include a ^10% systematic uncertainty in the absolute flux scale. 
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Table 2. Model Parameters 







fiducial 


high-g 


high-M* 


high-i 


hot 


non-Kep 


warp 


parameter 


units 
















log Niq 


[cm -2 ] 


19.00 


19.42 


19.28 


19.60 


19.00 


18.81 


18.83 


7 




0.99 


0.99 


0.66 


0.90 


0.99 


0.94 


0.93 


r c 


[AU] 


28 


28 


36 


24 


28 


32 


33 


T w (2-1) 


[K] 


77 


110 


75 


68 


77 


100 


100 


(3-2) 


[K] 


88 


115 


94 


90 


88 


99 


104 


q (2-1) 




0.38 


0.65 


0.39 


0.32 


0.38 


0.49 


0.53 


(3-2) 




0.44 


0.65 


0.51 


0.49 


0.44 


0.49 


0.53 




[m s _1 ] 


20 


20 


10 


10 


20 


20 


15 




[M ] 


0.8 


0.8 


1.5 


0.8 


0.8 


0.8 


0.8 


ho 


[°] 


5.8 


5.8 


6.0 


8.0 


5.8 


5.7 


7.5 


y 






















0.15 


n 


[AU] 












57 




X 



















0.15 





ST 




1 


1 


1 


1 


3 


1 


1 



Note. - The parameter values adopted in the modeling analysis in §3.2 and 3.3. 
Each column corresponds to a different model type, and each row represents a different 
model parameter (the subscript '10' denotes that parameter value at r = 10 AU). Note 
that only the "warp" model has a spatially varying disk inclination: in all other cases 
iio = i at all radii, and y = by definition (see §3.3.3). The parameter is only 
defined for the "non-Keplerian" model; in all other cases x = (or / = 1, at all radii; 
see §3.3.2). The parameter 5T corresponds to a constant scaling of the temperature 
profile for r < 4 AU in the "hot" model only: all other models have 8T = 1 by definition. 
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